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Abstract 

The multistep shell model was extended recently to incorporate both neutron and proton degrees 
of freedom and applied to study the structure N = Z systems with four, six and eight particles 



[arXiv:11 08.0269j. In this work we give a brief introduction to the formalism thus developed. A 
more detailed explanation with applications will be updated later. 

PACS numbers: 21.60.Cs, 21.10.-k, 21.30.Fe, 21.60.-n 
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I. INTRODUCTION 

In the multistep shell model method (MSM) ^ one solves the shell-model equation in 
several steps. In the first step one constructs the two-particle states. In the second step 
one proceed by solving the three- or four-particle states in terms of the two-particle states 
calculated in the first step. In our case we will solve the two-neutron plus two-proton system 
within a non-orthogonal overcomplete basis in terms of the (z^vr) (g) {un) excitations at the 
same time as the {uu) (tttt) ones. With the four-particle system thus evaluated, we will 
proceed to the next step and evaluate the six-particle system in terms of the four-particle 
states times the two-particle states. For the eight-particle system one can choose the MSM 
basis such that it consists of the four-particle states in the form (z/tt)® (z/tt) times themselves. 
Systems with more pairs can be described in the same fashion in successive steps. 

The Hamiltonian is given as 

i i 

+ ^{i3\ypn\kl)p\n]nipk + ^ ^{i3\Vpp\kl)p\p]piPk + ^ ^{ij\Vnn\kl)n\n]nink, (1) 

ijkl ijkl ijkl 

where pi and rti denote proton and neutron operators, respectively, e denote the single- 
particle energies, and V is the two-body interaction. The indices j, k, I label single- 
particle states. We will use the Greek letter 7„ to label the ra-particle np states. Since we 
will only consider cases with equal number of neutrons and protons outside a closed shell, n 
will be an even number such that the number of neutrons {n/2) is the same as the number 
of protons. Therefore the np states will be I72) = -P^(72)|0) where the np creation operator 
is -P^(72) = X^ij ^(^J! 72)^1^- III the same fashion the two-proton (two-neutron) creation 
operator is ^^(02) (-P^(/32))- The amplitude X of the two-body wave function is determined 
by solving the corresponding TDA equations, 

{W{a2) -Ep, - ep^){a2\p{pl\0) = ^^{PiPj\Vpp\piP2){a2\plp]\0), 

ViPj 
niUj 

(Vr(72) -Ep^ - en.^^l'i\p\n\\^) = '^{Pinj\Vpn\pin2){-f2\pln]\0). 

PiUj 

We have X = 1 for systems within a single-j shell. These correlated two-particle states are 
the basic building blocks in the construction of the MSM many-body basis vectors. 
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II. SYSTEM WITH TWO NEUTRONS AND TWO PROTONS 

The four-particle state is I74) = P''' (74)10) with 

Pt(^,) = Y,X{<^2^2;i,)P\a^)P\^2) + ^(727^; 74)^^(72)^^(72) (2) 

012P2 72<72 

where all possible like-particle and np pairs are taken into account. The physical meaningful 
quantities are the projections of the basis vectors upon the physical vector, which we denote 
as 

P(a2^2;74) = (74|P^(«2)P^(/92)|0), 

i^(7272;74) = (74|P^(72)P^(72)|0). (3) 
The orthonormality condition now reads 

5^,^. = X{a2/32; 74)^(^2/32; 7D + Y ^^2. 1^)FM2\ 74)- (4) 

Q2/32 72<72 

The norm of the MSM basis I7272) = (72)^^^(72)10), i-e., iV(7272;74) = vT^^^?^, may 
not be unity. 



II. 1. TDA Equation 

The dynamic matrix of the two-neutron two-proton system is given as 

(W^(74) - W^(72) - W^(72))(74|(P^(72)P^(72))74|0) = 
72^72" ^ PiP2nin2 ^2^2 



X(A1+A2) (74|(P^(72)^^(72"))74|0) 



^pi ^P2 ^ni ^n2 

) xB (74|(Pna2)P^(/32)),jO), 

02/32 V piP2nin2 ) 



and 

_ w{a2) - WiP2)){l4\iP\a2)P\fi2)hM 



Y{ E ^^^'^^ ^ ^^'^] ^ 7' " ~ ~ X ^ [(741(^^(7^0^-^7^)7.10). 



72^72" V PiP2nin2 '''2 '''2 
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where W denotes the corresponding n-particle energy. To obtain above equation we have 
assumed that the four-particle system was decomposed into two different blocks in terms of 

(tttt) (z^z^) and (z/vr) (z^tt). 

The A, B and C matrix elements are defined as, 



X X(pini;72)X(p2n2; 72)^(^1^2; 72)^(P2?ii; 72") ' 



Pi rii 72 

n2 P2 72 

72 72 74 



A„ — C_2^2pl+nl+n2+72+72"+74 



72727272" 



X X(pini;72)X(2)2n2; 72)^(^2^1; 72)^(^1^2; 72") < 



Pi ni 72 

^2 P2 72 
72 72 74 



Pi ni 72 



7272"2/32^(pini; 72)X(p2n2; 72)^(^1^2; a2)>^(nin2; /32) < 



P2 n2 72 

Q!2 132 74 



and 



Pi P2 Oi2 

€^ a2$2l2i2Y{piP2\a2)Y{nin2\l32)X{pini\i2)X{^^ ni n2 ^2 



I2 I2 74 



> . 



In all cases we use the same symbols to label states as well as the corresponding angular 
momenta. The coefficient Y is related to X by Y{ij\ a2) = (1 + dijY^'^Xiij] a2). 



II. 2. Overlap Matrix 

The overlap matrix is defined as follows, 

- (A1+A2) 

(0|(pt(^,)pt(^^))t^(pt(^^)pt(^^))^jO)= ® 

piP2nin2 

(0|(pt(«,)pt(^,));(pt(«^)pt(^^))^jO) 
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which correspond to the overlap between states of the forms {vn ® vt:\vt: ® utt), {vn ® 
vt:\vv (g) tttt) and {vv ® t:t:\vv (g) tttt), respectively. 

III. SYSTEM WITH THREE NEUTRONS AND THREE PROTONS 

For the six-particle case we will use the MSM partition of two- times four-particles. Thus 
the corresponding wave function will be |76) = P^(76)|0), where 

P\^,) = J] X(7274; 76)^^(72)^^(74). (5) 

7274 

and 

|76) = Xl"^*^^274;76)(7274|76) X]^(a2^2;74)(tt2/32|74) X]^(Pini;72)(Pini|72) 



7274 a2P2 Pin-l 



JX]^(P2]?3;a2)(p2P3|a2) X \ ^Y{n2n2]l32){n2nz\(52) Pin[p\pln\,n\\Q). 



2^ V-r^x-o, z./vr^x'oi ^1 2 
P2P3 n2n3 

As before, we will evaluate the projection of the basis vectors upon the physical vectors, i.e., 
F(7274; 75). In this six-particle case one can also view the MSM basis elements as the direct 
tensorial product of three pairs which takes the forms vt: (S) uti ® i^tt and vt: ® ® tttt. 

III.l. TDA Equation of the 2 x 4 Block 

For the partition of one np pair times the 4-particle system, the dynamic matrix is given 

as 

(W^(76) - W{^2) - W^(74))(76|(7l7l)76|0) = 
E{ E E {W{i2) + W{P',)-e,,-e^,-er,,-e^,)xA, 

7274 ^ Pi"i"2n3 a2/32/3^6'4 

>l (6) 

+ E E (^(^2) + W(a'^) - Sp, - Sp, - - X A2 |(76|(72S;^)76|0), 

piP2P3ni a2/92a204 

where 

Ai = (-l)"^+"3+^^+^2+^2+74+74^2^2;^^^/^^^/^2 

xX{pini;j2)y{n2ns; /32)X{pin3; j2)y{nin2; /3'2)X{a2/32; j4)F{a2/32, 74) 



Pi ni 72 



X < 



ns n2 P2 
72 P2 0, 



72 (32 0, \ \ 7^ (3', 9, 



«2 76 74 "2 76 74 



5 



A2 = (-l)"l+^^+^^+"^+^^+'^^72«272«'27474^4 

xX(pini; 72)1^ (P2P3; a2)X{p3ni; 72)^(^1^2; a2)-^("2/32; 74)i^(a2/52; 74) 

Pi rii 72 

X < p2 P3 ^2 
a'2 72 04 

The overlap matrix of 2 x 4 block is given as 



72 Q;2 04 M 72 «2 04 
/52 76 74 I 1 /52 76 74 



+ E E E E 

p\nin2nz 02/32/32^4 piP2P3ni a2/32a2</'4 

+ E E » 

PiP2P3nin2n3 a2/32a2/32V'2 



, where 



B = (-l)^^+^^+^H2«2/3272«2/32747;^2 

xX(pini; 72)F(p2P3; Q;2)l^(n2n3; P2)X{psn3; 72)5^(^1^2; agjl^Kna; /^a) 

72 V^2 74 

xX(q;2^2; 74)^(0:2^2; 74) ^ P2 ^^2 -02 Pa JT-a 72 !" ^ 

72 76 74 







> 










1" 


ni 


72 




V2 


^2 


-02 


P2 


n2 


-02 


> < 


Pa 


JT-a 


7^ > 


«2 

V 


^'2 


7^ 




a;2 

V. 


^2 


74 



(8) 



(9) 



(10) 



III. 2. Transformation to the 2x2x2 block 

The transformation from the 2x4 block to the 2x2x2 block is given as 

(76|(7l«^/3l)76|0) = E(^6|(7l7l)76|0)(74|(4/32^)74|0), 

74 

(76|(7h2'7f )76|0) = E(^6|(7hl)76|0)(74|(7^^7f )74|0). 

74 

The overlap for the 2x2x2 coupling is 
{^\{l\AB\%{i2'c\D\\,\^^ 



747i 



IV. FOUR-PROTON FOUR-NEUTRON SYSTEM 



We will describe the eight-particle states as I78) — P''^(78)|0), where 



^^(78)= ^(«4/34;78)P^(«4)P^(/94). 
a4</94 



(11) 
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and 



Its) =E^(^474;78)(7474l78) E^(«2^2;74)(«2^2|74) E^(«2^2;7D(«2^2l7^ 
74t4 azfe 



piP2 nin2 



X -^Y{p3Pi;a2){p3P4\a'2) x -^Yinsn^; f32){n3n4\P2)p\pWinlplplnln\\0) 



P3P4 n3T»4 



The dynamic matrix for the 4x4 block is defined as 



(M^(78) - W^(74) - W^(7D)(78|(747D7s|0) = E Y^TfV^, ^ 



^(74) + ^^(74") - ^(«2) - W{^i) - Wio!^) - j X (Ai + A2) 

+ 5Z 5Z (^("2) + ^(«2') - %i - - - ^4) (Kl + K2) 

+ E E (w^(/320 + ^(/320-£ni-£n2-£n3-£n4)(Cl + C2)l(78|(7:^^^^^^ 



where 



A, = (_l)/32+/3^+74+7r X 74747474" 



^2 /32 74 

xX(«2/32;74)^(«2/32;7;)i^(a2/32;7:)i^(«2/32;74")^ 74 



74 74 78 



A2 = (-1) 



/32+/34+74+74"+78 



X 747^7474 



X x(a2/32; 74)^(^2/32; 74)^(«2/32;74')i^(«2/32; 74') { 



«2 f^2 74 

/S^ 0^2 7i 



74 74 78 
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Bi = J](-l)a2«^«'2'«r74747474"^l04 

X X(a2^2; 74)^(«2^2; 74)^(«2^2; ^'DF{a'^P',; ^T)Y{prp,; a2)Y{pm; ^2) 





f 


Pi 


P2 


OL2 


> 




Ci2 


/32 


74 






/32 


a ^ 


X r (piPs; ai2)^(P2P4; ^2) < 




P3 


Pi 


"2 




> < 


0^2 




7i 


> < 


< 




7r ^ 






a'i 


«r 




> 




^4 


04 


78 

J 




^4 


04 


78 

J 



®2 = X a2a'2ala'm'AlllT0Ul 

044)4 

X X(a2^2; 74)^(«2^2; 74)^(«2^2; ll)F{a'^P2; l7)Y{piP2; a2)Y{p3P4; «2) 





f 


Pi 


P2 


Ci2 


> 




Ci2 


P2 


> 

74 






(^'2 




X r (pips; «^')^(P2P4; aa") < 




P3 


P4 


«2 




> < 




P'2 


7i 




a- 


P2 


7r 










^4 


> 




< 


04 


78 

J 




. ^4 


04 


78 J 



Ci = 5^(-l)^2^^^^'^r747;7:74'^l0l 

X X{a2f32; 74)^(0;^^^; 7D^(«2^2 ; 11)^0^2^2; 74")^(^i^2; ^2)1^(^3/14; ^2) 





Til 


n2 


/92 ^ 




Cl2 


> 

/32 74 




OL2 




7;' ^ 


xF(nin3;/3^')>^Kn4;/3^")< 








> < 


«2 


/32 74 


> < 






71" > 












^4 


04 78 

) 




< 


04 


78 



C2 = E(-l)'''"^"^^''^2^2«'747;7:7r^l04 

X X(a2^2; 74)^(«^^^; i:)F{ol^&h il)F{a2^'^; l7)yinin2; ^2)1^(^3/14; ^^) 





rii 


n2 


/32 ^ 




Ci2 


/32 


> 

74 






/32' 


7l' ^ 


xY{nins;/3!^)Y{n2n^;/3!^')< 


ns 




^2 


> < 


«2 




7i 


> < 


OL2 




iH > 








) 




^4 


04 


78 

) 




^4 


04 


78 
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IV. 1. Overlap Matrix 



The overlap matrix of 4 x 4 block is given by 

+ E E E (®i+®2)+ E E E c^i+'c^) 

02/3202/32 OjOj P1P2P3P4 02/3202/^2 /32'/32" "I"2n3n4 

(-'-■^J 

+ E E EE". 

02/3202/32 ai^^i^ai^'^!^' P1P2P3P4 nin2n3n4 



where 



4^4 



D = 5]d2/92«2/52«2/52«2'/52'747i747r^: 

^4(^14 

XX(a2;52; 74)^(«2/32; 7D^(«2/32 ; 7:)^(«2"/32"; 7r)>^(PlP2; «2)1^(P3P4; «^ 

xF(piP3; a2)>'(P2P4; a2")>^(^i^2; /32)V'(n3n4; f3',)Y {n.ns, P'^)Y{n,n,; f3'^' 



X < 



Pi P2 0L2 
PZ Pi "2 



> < 



rii 712 /S2 
na Hi 
P'i 04 



> < 



"2 /S2 74 
«2 ^2 74 
6*4 04 78 



> < 



a'i a 

a- 

^4 04 78 



> . 



(13) 



IV.2. Transformation to the 2x2x2x2 coupling 



The transformation from the 4x4 block to the 2x2x2x2 block is given as, 



^741 



(78|(«2/S2a^/S2)78|0) = $^(78|(747;)78|0)(74|(a2/S2),4|0)(7il(«2/92)' 

747^ 

(78| («2/327272)78 10) = 5^(781 (747^78 10) (74| (^2/32)^4 10) (74I (7272)74 10)' 

7474 

(78|(72727272078|0) = $^(78|(747D78|0)(74|(7272)74|0)(74l(^^^^^^^ 



7474 



The overlap matrix is as follows, 



(14) 



(0|(A2S2C2L'2);(£;2F2G'2i/2),s|0) =E E(74|(^2S2),J0)(7;|(C2L'2),4|0) 



7474 74'74" 



x(0|(747;);(7M78|0)(7^1(£;2F2),4'|0)(7ri(G2//2),4"|0). (15) 
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V. SUMMARY 



Summarizing, in this work we introduced the formahsm to study the structure of the 
wave function ol N — Z systems with two (Section II), three (Section III) and four (Section 
IV) np pairs. Systems with more pairs can be described in the same fashion in successive 
steps. 
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